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Abstract 

Supersymmetric grand unified models based on the gauge group SU(5) often require in ad- 
dition to gauge coupling unification, the unification of 6-quark and r-lepton Yukawa couplings. 
We examine SU(5) SUSY GUT parameter space under the condition of b — r Yukawa coupling 
unification using 2-loop MSSM RGEs including full 1-loop threshold effects. The Yukawa-unified 
solutions break down into two classes. Solutions with low tan/3 ~ 3 — 11 are characterized by 
nig ~ 1 — 4 TeV and mq ~ 1 — 5 TeV. Many of these solutions would be beyond LHC reach, 
although they contain a light Higgs scalar with nih < 123 GeV and so may be excluded should 
the LHC Higgs hint persist. The second class of solutions occurs at large tan /3 ~ 35 — 60, and are 
a subset of t — b — r unified solutions. Constraining only b — r unification to ~ 5% favors a rather 
light gluino with rrig ~ 0.5 — 2 TeV, which should ultimately be accessible to LHC searches. While 
our b — t unified solutions can be consistent with a picture of neutralino-only cold dark matter, 
invoking additional moduli or Peccei-Quinn superfields can allow for all of our Yukawa-unified 
solutions to be consistent with the measured dark matter abundance. 
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1 Introduction 

Grand unified theories (GUTs) based upon the Lie group SU(5) are very compelling in that they 
unify the disparate gauge groups of the Standard Model (SM) into a theory based upon just a 
single gauge group [1]. Furthermore, the SU(5) theory provides a rationale for the seemingly ad- 
hoc weak-hypercharge assignments of SM matter fields. A third triumph occurs in that models 
based on SU(5) unification unify the 6-quark and r-lepton Yukawa couplings at the unification 
scale Mqut ~ 10 16 GeV: then, renormalization group effects provide roughly the correct values 
of nib and m r at low energy scales [2] . 

Adding supersymmetry to SU(5) grand unified theories seems essential in order to stabilize 
the vast hierarchy separating the weak scale from the GUT scale [3]. The additional SUSY 
degrees of freedom contained in the Minimal Supersymmetric Standard Model (MSSM) alter the 
RG running of gauge couplings below the scale Mqut- Indeed, the celebrated unification of the 
three SM gauge couplings within the MSSM is often touted as indirect evidence for weak scale 
SUSY, and for the possible existence of a SUSY GUT theory [4] . 

Some drawbacks to the SUSY SU (5) theory include the incorporation of a rather awkward set 
of Higgs multiplets which are necessary for appropriate GUT symmetry breaking. Foremost among 
these problems is embedding the MSSM Higgs doublets H u and Hd into a 5 and 5* respectively 
of SU(5): one must then explain why the color triplets obtain GUT scale masses while the MSSM 
multiplets receive weak scale masses: the so-called doublet-triplet splitting problem. In addition, 
in SUSY GUT theories based in four spacetime dimensions with spontaneous GUT symmetry 
breaking via the Higgs mechanism, protons are expected to decay even in the SUSY theories with 
rates which now seem excluded by experiment [5, 6]. 

However, if one formulates SU(5) SUSY GUT models in five or more spacetime dimensions, 
then the GUT symmetry can alternatively be broken by compactification of the extra dimensions 
on an appropriate manifold such as an orbifold [7]. This method of symmetry breaking can 
solve the doublet-triplet splitting problem while suppressing or even eliminating proton decay. 



Construction of SUSY GUT models in extra dimensions can solve many problems endemic to 
4-d theories while preserving many of the compelling features of unified theories. Indeed, in 
string theory it is possible for SUSY GUT theories to emerge on four or more dimensions as the 
low energy (GUT scale) effective theory, where the 6-7 additional stringy dimensions must be 
dispensed with anyway [8]. 

In this paper, we seek to avoid the very model-dependent physics associated with the GUT 
sector and possible extra dimensions, and use instead data and some general SU(5) SUSY GUT 
characteristics as a guide to what weak scale physics should look like at colliders such as the 
LHC We will assume here that nature is described by an SU(5) SUSY GUT theory at energy 
scales Q ~ M GUT ~ 2 x 10 16 GeV. At Mgut, the MSSM superfields Q, U c and E c live in the 
antisymmetric 10 of SU(5): ip 13 , while the D c and L superfields live in the 5*: fa. (Here, i and 
j are SU(5) indices running from 1-5.) The MSSM Higgs doublet is an element of a 5*: Hu, 
while H u is an element of a 5: H^. The superpotential / then contains the terms [9, 10] 

/ 3 \fttijklmftH kl Hl? + V2f b ^faH Xj + ii H H u H* + ■■■ (1) 

so that the third generation Yukawa couplings fb and f T are unified at Mqut, but are distinct 
from f t . 

The soft SUSY breaking terms in an SU(5) SUSY GUT theory are expected to include: 
£ so ft 3 -m 2 Hl \Hi\ 2 - m 2 H2 \H 2 \ 2 - ml\<f>\ 2 - ml tr{ip^ip} - ^m 1 / 2 X a Xa 



\A t f t e ijklm ^i> kl H? + V2A b f h ^faH l3 + h.c. 



(2) 



For SU (5) SUSY GUT models with b — r Yukawa coupling unification, we will adopt a GUT 
scale parameter space given by 

m 5 , raio, , rn 2 Hu , rn 2 Hd , m 1/2 , A t , A b , tan/3, sign(fi) (3) 

where we identify m 2 H ^ = m 2 H2 and m^- = m 2 Hi . We also take the top quark pole mass to be 
nit = 173.3 GeV, in accord with recent measurements from CDF and DO [11]. 

Recent previous work on Yukawa coupling unification has focused on t — b — r unification 
which is expected to occur in the simplest 50(10) SUSY GUT models [12]. In these models, 
unification of all matter (super)fields of a single generation into a 16-dimensional spinor tp occurs. 
The two MSSM Higgs multiplets are also unified into a 10-dimensional Higgs representation <fi. In 
these models, it was found that t — b — r Yukawa unification can occur if the soft SUSY breaking 
(SSB) parameters are related as A^ ~ 2mf ~ 4mf 6 [13, 14, 15, 16]. With matter scalar SSB 
masses in the multi-TeV range and gaugino mass Wjm as small as possible, these relations lead 
to a weak scale sparticle mass spectrum of the inverted mass hierarchy type [19]: first/second 
generation squarks in the 5 — 20 TeV range while third generation scalars are at < 1 TeV, as 
required by naturalness. Either a "just-so" splitting of Higgs SSB terms (with rrijj > rn 2 Hu at 
the GUT scale), or D-term splitting of scalars (in the DR3 model [20]) is required for radiative 
electroweak symmetry breaking (REWSB) [21]. 

In SUSY models with t — b — r Yukawa unification and unified gaugino masses, there is a 
tendency in the sparticle mass spectrum for rather light gluinos with nig < 500 GeV [13, 15, 16] 
(although solutions can also be found with significantly higher gluino masses [22]). Recent searches 
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for gluino pair production in Yukawa-unified models require rrig > 500 GeV, placing some stress 
on this class of models [23]. One path to relieve such stress is to assume a two-stage breaking 
pattern 

50(10) -»• SU(5) -»• SU(3) C x SU(2) L x U(1)y- (4) 

In this case, one might expect a high degree of b— r Yukawa unification, but perhaps a lesser degree 
of t — b — t unification. If the scales of the two stages are significantly separated, then one needs to 
take into account the evolution of the SSB terms above Mqut- Such super-GUT effects can lead 
to sufficiently different sparticle spectra with interesting phenomenology [9, 10, 17]: for example, 
the no-scale scenario can be made compatible with experimental constraints [18]. The breaking 
pattern (4) may also allow for heavier gluino masses to occur in the range of m g ~ 0.5 — 1 TeV. 
Such gluino masses should be accessible to LHC SUSY searches with ^/s = 7 TeV and 20 — 30 fb _1 
of integrated luminosity [24] . 

Spurred by these developments, the authors of [25, 26] investigated t — b — r Yukawa unifi- 
cation in the framework of SUSY SU(A) C x SU(2) L x SU(2) R [27] (4-2-2, for short). The 4-2-2 
structure allows one to have non-universal gaugino masses while preserving Yukawa unification. 
An important conclusion reached in Ref. [25] is that with the same sign but non-universal soft 
gaugino masses, Yukawa unification in 4-2-2 for ii > is compatible with neutralino dark matter, 
with gluino co-annihilation [25, 28] playing an important role. 

By considering opposite sign gauginos with fj, < 0, Mi < 0, M3 > 0, (where M2 and M3 are 
the SSB gaugino mass terms corresponding to SU{2)i and SU(3) C respectively, it is shown in 
Ref. [26] that Yukawa coupling unification consistent with the experimental constraints can be 
implemented in 4-2-2. With < and opposite sign gauginos (M2 < 0, M3 > 0), Yukawa 
coupling unification is achieved for rriiQ > 300 GeV, as opposed to mi6 > 8 TeV for the case of 
same sign gauginos. The finite corrections to the b-quark mass play an important role here [26]. 
Note that with M2 < 0, M3 > and \i < 0, we can obtain the desired contribution to (g — 2)^ [29]. 
This enables us to simultaneously satisfy the requirements of t — b — r Yukawa unification in 4-2-2, 
neutralino dark matter and {g — 2)^, as well as a variety of several other bounds. 

Encouraged by the abundance of solutions and co-annihilation channels available in the case 
of Yukawa unified SUSY 4-2-2, Yukawa unification in SO(10) GUT was explored in [30] with 
non- universal MSSM gaugino masses at Mqut- This scenario can arise from non-singlet F-terms, 
compatible with the underlying GUT symmetry [30]. Furthermore, the soft masses for the two 
scalar Higgs doublets are set equal {mH n =fnH d ) at Mqut- It is intriguing to note that in these 
models, rather precise t — b — t Yukawa unification also happens to yield a mass for the lightest 
CP-even Higgs boson in the 122 — 124 GeV range [13, 31]. There is an approximately 2 GeV 
theoretical uncertainty in this calculation. 

The remainder of this paper is organized as follows. In Sec. 2, we outline details of our 
sparticle mass spectra calculation, along with the requirement of b — r Yukawa unification, and 
constraints from earlier collider and B decay searches. In Sec. 3, we show preferred SU(5) model 
parameter choices which lead to b — r Yukawa unification. The solutions divide into two classes: 
1. those with low tan/3 ~ 3 — 11 for which the top Yukawa coupling /(>/(, ~ f T at the GUT 
scale, and 2. those at high tan/3 ~ 35 — 60, which give t — b — r Yukawa quasi-unification, i.e. 
ft ~ fb — fr a t Q = Mqut- The low tan/3 solutions may be eliminated if the LHC hint of Higgs 
at irih — 125 holds true. Otherwise, LHC direct searches for sparticles can only cover a portion 
of the low tan f3 solutions, while searches for gluino pair production at LHC can cover nearly all 
of the high tan /3 solutions. In Sec. 4, we discuss aspects of the relic abundance of dark matter 
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for b — t Yukawa-unified models. While neutralino-only dark matter can be accommodated by a 
variety of co-annihilation, resonance annihilation or higgsino annihilation processes, generically we 
expect a standard overabundance of neutralinos. Either an overabundance or an under abundance 
of neutralino dark matter can be brought into accord with the measured CDM abundance by 
invoking either additional late-decaying scalar (moduli) fields, or by invoking a Peccei-Quinn 
axion superfield (containing axion, saxion and axino components) which is needed anyway as a 
solution to the strong CP problem. A summary and conclusions are presented in Sec. 5. 

2 Calculation of sparticle mass spectra with b — r Yukawa unifi- 



For our calculations, we adopt the Isajet 7.80 [32, 33] SUSY spectrum generator Isasugra. Isasugra 
begins the calculation of the sparticle mass spectrum with input DR gauge couplings and /&, f T 
Yukawa couplings at the scale Q = Mz {ft running begins at Q = mt) and evolves the 6 couplings 
up in energy to scale Q = Mqut (defined as the value Q where g\ = #2) using two- loop RGEs. 
We do not strictly enforce the unification condition 53 = 51 = g2 at Mqut, since a few percent 
deviation from unification can be assigned to unknown GUT-scale threshold corrections [34]. At 
Q = Mqut, the SSB boundary conditions are input, and the set of 26 coupled two-loop MSSM 
RGEs [35] are evolved back down in scale to Q = Mz- Full two- loop MSSM RGEs are used for 
soft term evolution, while the gauge and Yukawa coupling evolution includes threshold effects in 
the one-loop beta-functions, so the gauge and Yukawa couplings transition smoothly from the 
MSSM to SM effective theories as different mass thresholds are passed. In Isajet 7.80, the values 
of SSB terms which mix are frozen out at the scale Q = Msusy = yJ m t L m t R i while non-mixing 
SSB terms are frozen out at their own mass scale [33]. The scalar potential is minimized using 
the RG-improved one-loop MSSM effective potential evaluated at an optimized scale Q = Msusy 
which accounts for leading two- loop effects [36]. Once the tree- level sparticle mass spectrum 
is computed, full one-loop radiative corrections are calculated for all sparticle and Higgs boson 
masses, including complete one- loop weak scale threshold corrections for the top, bottom and 
tau masses at scale Q = Msusy [37]. These fermion self-energy terms are critical to evaluating 
whether or not Yukawa couplings do indeed unify [38]. Since the GUT scale Yukawa couplings 
are modified by the threshold corrections, the Isajet RGE solution must be imposed iteratively 
with successive up-down running until a convergent sparticle mass solution is found. For most 
of parameter space, there is excellent agreement between Isajet and the SoftSUSY, SuSpect and 
Spheno codes, 6 although at the edges of parameter space agreement between the four codes 
typically diminishes [39]. 

We searched for Yukawa-unified solutions in the SU(5) parameter space (3) in two stages. 
First, we performed the MCMC scan [16, 40] over the large parameter range 



6 These three codes invoke an "all-at-once" transition from MSSM to SM effective theories in contrast to the 
fsasugra approach. 



cation 



mio, rn 5 , m Hu , m Hd : - 20 TeV, 
m 1/2 : - 2 TeV, 
-60 TeV < A t , A b < 60 TeV, 
tan/3 : 1.1 - 60. 



(5) 
(6) 

(7) 
(8) 
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We identify several solutions with good b — r Yukawa unification [R < 10%) and the neutralino 
Relic Density ^o/i 2 within the WMAP bound [41]. Those solutions become centers of second- 
stage scans with narrower parameter ranges, where we look for more solutions with good b — r 
Yukawa unification, good ^^oh 2 and try to make the spectra as light as possible. Each generated 
parameter set is entered into Isasugra using the non-universal SUGRA model inputs, and our 
initial selection criteria is that the points generate a neutralino \\ as lightest MSSM particle, and 
appropriate REWSB. In plots to follow, these points are labeled with gray color. 

We next require the following bounds (inspired by LEP2/Tevatron searches) on sparticle 
masses: 

• m h > 114.4 GeV, 

• mi v mg i > 100 GeV, 

• m fl > 105 GeV, 

• m-g > 250 GeV, 

• m-± > 103 GeV. 

Xi 

Using Isatools [42, 43] and [44], we also require the following bounds from heavy flavor (i?-physics) 
to be respected: 

• BF(B S -> n+fi-) < 1.1 x 10~ 8 [45], 

• 2.85 x 10" 4 < BF(b -> s>y) < 4.24 x 10~ 4 [46], 

• 0.15 < BF SUSY (B U -> tu t )/BF sm (B u tv t ) < 2.41 [46]. 

Points passing both mass and -B-physics cuts are labeled as red or green. 

For each solution, we calculate the degree of b — r Yukawa unification at the GUT scale via 
the i?-parameter: 

max{f b , f T ) 
min(f b ,f T ) 

Thus, R = 1.0 would tag a solution with perfect b — r Yukawa coupling unification. 



3 Results 

3.1 SU(5) parameters required by b — t unification 

Our first results are shown in Fig. 1. In Fig. la), we plot the value of R vs. tan/3. Points with 
R < 1.05 have a high degree of /& — f T Yukawa unification. We see immediately that the solutions 
break up into two classes: low tan (3 ~ 3 — 11 and high tan (3 ~ 35 — 60. We color code the resulting 
points according to this criteria: red points have tan/3 < 20 and green points have tan/3 > 20. 
Apparently no b — r unified points can be generated for 11 < tan/3 < 35 which satisfy the mass 
and S-physics bounds. 

We will first discuss the low tan/3 ~ 3 — 11 solutions. From Fig. lb) and c), we see that this 
class of solutions requires mio ~ 1 — 4 TeV and ~ 1 — 5 TeV. In addition, from frame d), 
we see that rather large values of ~ 0.6 — 1.5 TeV are required. Such high m^ 2 values 
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Figure 1: b — r Yukawa unificaion as function of model parameters for [i > and m% = 173.3 GeV. 
Gray points satisfy REWSB and neutralino as LSP conditions. Red and green points satisfy 
additional spartcile mass and S-physics bounds and have tan j3 < 20 and tan/3 > 20, respectively. 
The horizontal dashed line indicates the 5% Yukawa unification. 
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lead to gluinos with m g > 1.3 TeV. Moreover, as noted above, if we insist that R < 1.05, there 
is a clear distinction between low and high tan /3 solutions. However, if we relax this condition 
(R > 1-3), this distinction between the low and high tan/3 solutions disappears. Another point 
to note is that here we have b — r unified solutions for low tan/3. Such solutions are ruled out 
in the CMSSM because of the Higgs mass bound [47]. In our scans we find solutions with tan/3 
as low as 3.5. This limit can be slightly changed by varying mt and nib- From frames e) and 
f), we see that R < 1.05 solutions require the GUT-scale trilinear SSB parameters to be in the 
range — 2.3mio < A t < and Af, ~ (10 — 20)ms. The latter parameter A^Mqut) again strongly 
differentiates between the low and high tan /3 solutions. The low tan (3 solutions are characterized 
by fb — fr ~ 0.04 at the GUT scale, while ft(McuT) ~ 0.5, i.e. there is a large disparity between 
ft and fb ~ f T . This class of models might be indicative of an SU(5) GUT theory, but with no 
connection to 5O(10) (unless we proceed to 5O(10) models where the MSSM Higgs doublets live 
in separate 10s of 5O(10) [10]). 

In contrast, from Fig. 1 we see the green points with tan /3 ~ 35 — 60 require mio ~ m§ ~ 
5 — 20 TeV, i.e. these solutions require multi-TeV matter scalars just as do SUSY models with 
t — b — t Yukawa unification [15, 16, 48]. Furthermore, we see from frame d) that the high tan/3 
points require much lower values of m-1/2 ;$ 0.7 TeV, which leads to an upper bound on the gluino 
mass of m g < 2 TeV. From frame e), we find that — 2.8mio < At < — 1.8mio, with also a few 
solutions around At ~ 2.2mio- In frame /), we see that Ab is much less correlated, with \ A^\ < 3ms. 
The high tan/3 ~ 45 — 55 solutions also tend to have a high degree of ft — fb — fr ~ 0.55 
unification, whereas the solutions with tan /3 ~ 35 — 45 tend to have ft ~ 0.55, but fb ~ f T ~ 0.47. 
These latter solutions might be indicative of an 5O(10) SUSY GUT which has broken to SU(5) 
at a higher mass scale than where SU(5) breaks to the SM gauge group. 

In Table 1, we list low tan/3 and high tan/3 benchmark solutions for illustration. Points 1 
and 2 belong to the set of low tan /3 Yukawa unified points and also represent ^4-resonance and 
sbottom co-annihilation scenarios. Points 3 and 4 are representative of high tan /3 Yukawa unified 
solutions, where point 3 depicts a stop co-annihilation solution, while point 4 shows a large O/i 2 
value. 

3.2 SU(5) preferred masses 

We next proceed to examine some derived parameters associated with the Higgs/higgsino sector 
from SUSY models with b — r unification. In Fig. 2 a), we show the correlation of the degree of 
the Yukawa unification with the parameter [i. The magnitude of [i is determined by minimization 
conditions on the Higgs scalar potential. Here, we see that the low tan /3 solutions also give rise 
to a range [i ~ — 5 TeV: i.e. \i is bounded from above, and furthermore \x can be well below 
the 1 TeV scale. This may allow for the lightest neutralino x? to be of mixed bino-higgsino type, 
which gives rise to WMAP-allowed values of thermal neutralino relic abundance. However, the 
high tan /3 solutions require rather large values of fi, typically in the multi-TeV range, so that for 
these solutions we would expect the Xi state to be a nearly pure bino. Here too we can see the 
separation of low and high tan/3 solutions for R < 1.05. 

In frame b), we show correlations of the /^-parameter with the mass of the CP-odd Higgs 
boson A. The low tan /3 solutions require tjia ~ — 6 TeV. The rather low range of rriA may 
allow for neutralino annihilation through the ^4-resonance in the early universe. Note that for low 
tan /3 values, we can have LHC accessible solutions for txia if we require b — r unification better 
than 5%. Meanwhile, the high tan/3 solutions tend to have rriA inhabiting the multi-TeV range, 
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Figure 2: Rb T versus parameters in the Higgs sector for /i > and nit = 173.3 GeV. The color 
coding is the same as in Fig. 1. 
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Table 1: Input parameters and resulting mass spectra and rates for several sample points from 
the scan. All masses and dimensionful parameters are in GeV units. 

so that A resonance annihilation is unlikely to be a possibility. 

In frames c) and d), we plot R as function of the mass of the light CP-even Higgs boson h 
for low tan /3 (red) and high tan f3 (green) solutions, respectively. In this case, we see that the 
low tan/3 solutions require < 123 GeV, with usually much lower. Meanwhile, the large 
tan ft solutions allow for rrih ~ 123 — 133 GeV. Recently, some evidence has been reported from 
the Atlas and CMS experiments [49, 50, 51] for a SM-like Higgs boson very near to ~ 125 GeV. 
If this result is maintained by the factor of 4-6 more data from LHC expected in 2012, then it 
would likely rule out the low tan ft b — r unified solutions, while maintaining consistency with the 
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high tan /3 solutions! 

Because good b — r Yukawa unification requires large values of SSB masses, 771,5 ^ 1 TeV, first 
and second generation sleptons are rather heavy with masses greater than ~ 1 TeV (> 4 TeV for 
high tan/3 values). This leads to the large suppression of the SUSY contribution to the muon 
anomalous magnetic moment that arizes at 1-loop level from diagrams involing smuon and 

muon sneutrino. As result a^ USY values, that we computed using the IsaAMU [52] subroutine, is 
always several orders of magnitude below the extracted discrepancy Aa^ = (28.7±8.0) x 10 -10 [53] . 
This can be seen in several sample points we listed in Table 1. 

3.3 Prospects for LHC SUSY searches 

In this section, we discuss prospects for detection of b — r Yukawa-unified SUSY at the CERN 
LHC pp collider with either \fs = 7 or 14 TeV. In Fig. 3, we show solutions which pass the mass 
and -B-physics cuts - but also with R < 1.05 - in the vs. m g plane. The value of is 
meant to exhibit a typical first/second generation squark mass. The red low tan/3 points all have 
777^ > 1 TeV with rriq ~ 1.5 — 5 TeV. In this case, we find rriq ~ rrig or slightly heavier. The reach 
of LHC with yfs = 7 TeV and 20 fb" 1 (LHC7) extends to m g ~ rriq ~ 1.5 TeV [24], while LHC 
with y/s = 14 TeV and 100 fb" 1 extends to nig ~ rriq ~ 3 TeV [54]. Thus, about half the red 
points will be within reach of LHC14, while those with rriq > 3 TeV and rrig > 2 TeV will likely 
be beyond LHC14 reach. LHC luminosity or energy upgrades will be necessary to probe more 
deeply into the low tan /3 Yukawa-unified space of solutions. 

In the case of high tan/3 solutions (green points), we see that rriq exceeds - and frequently 
far exceeds - 5 TeV. Meanwhile, the gluino mass is bounded from above, with rrig almost always 
< 2 TeV. In this case, LHC searches will focus on gluino pair production [55]. For rriq 3> m g , 
the LHC7 reach is to rrig ~ 1 TeV [24], while LHC 14 reach is to rrig ~ 1.7 TeV [54]. Recent 
work on LHC signatures in the case where rriq 3> m g implies the maximal LHC14 reach in models 
with gaugino mass unification is to 777^ ~ 2 TeV in the XF~X2 ~~ ^ + Ey? 88 channel [56]. Thus, 
we expect LHC7 to probe the region rrig < 1 TeV, and LHC 14 to nearly cover the remaining 
parameter space with 100-1000 fb _1 of integrated luminosity. 

4 Dark matter relic density in b — r unified models 

It has been noted long ago that t — b — T Yukawa-unified models tended to give a huge overabun- 
dance of thermally produced neutralino-only dark matter [15]. 7 To this end, we adopt the IsaReD 
relic density calculator [59] to compute the thermally produced neutralino abundance ^o^ 2 from 
b — t unified models. The value of f^o/7 2 versus m^o is plotted for solutions with R < 1.05 in 
Fig. 4. The red solutions with low tan/3 span a range 10~ 3 < f^ofr 2 < 10 2 . The large tan/3 
solutions populate a much larger range of 10 -3 < f^o/7 2 < 10 5 owing to the suppression of y? 
pair annihilation by exchange of multi-TeV scalars in the relevant Feynman diagrams. This is to 
be compared with the CDM abundance tt C DMh 2 = 0.1109 ± 0.0056 reported by WMAP7 [41], 
which we indicate as a horizontal black line. 

7 In Ref. [57] , it is suggested that A-resonance annihilation may be available for bring neutralino-only CDM into 
its measured range. However, a number of other authors have failed to reproduce t — b — T unified solutions with 
very low uia values such that 2m^o ~ jua [15, 16, 58]. 
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Figure 3: Distribution of points from the scan in the mass plane of gluino and lst/2nd generation 
squarks for [i > and mt = 173.3 GeV. Gray points satisfy REWSB and neutralino as LSP con- 
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and represent tan/3 < 20 and tan/3 > 20, respectively. The dashed lines represent approximate 
reaches for LHC7. 
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Figure 4: Neutralino relic density fi^o/i 2 versus the neutralino mass m^o from the scan with 
fi > 0. All points satisfy mass bounds, -B-physics bounds and have R < 1.05. Red and green 
color represent solutions with tan/3 < 20 and tan/3 > 20, respectively. The solid horizontal line 
represents the WMAP measured value [41]. 
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For the low tan f3 case, it is possible to gain solutions with ^4-resonance annihilation, higgsino 
annihilation or stop, sbottom or stau co-annihilation. For these processes to be significant, certain 
mass conditions needs to be fulfilled: the mass gap between Xi an d the next ligtest sparticle needs 
to be within ~ 15% for coannihilation or m^o — 2mA for A-resonance annihilation. For higgsino 
annihilation, x± needs a sizable higgsino content, which makes it close in mass with the chargino 
Xi which is a higgsino-wino mixture. This is illustrated in Fig. 5, where we show solutions 
with R < 1.05 in the a) vs. m^o plane, the b) vs. m^o plane, c) the vs. m^o 
plane, d) the m~± vs. m^o plane and e) the wa vs. m^o plane. In each case, the approximate 
needed conditions for co-annihilation, resonance annihilation or mixed higgsino annihilation are 
illustrated by diagonal black lines. The neutralino-sbottom co-annihilation solutions shown in 
Fig. 5 c) are consistent with the results presented in Ref. [60], where it is shown that it is not 
trivial to have such a scenario. One needs to go to the more flexible SU(5) GUT scale boundary 
conditions to realize such solutions. 

For the high tan/3 case, the value of Q^oh 2 tends to be very high, well over fl^oh 2 ~ 1. We 
see that there are no solutions with sbottom and stau coannihilation, as shown in Figs. 56) and 
c). This is because Yukawa unification requires Af, = A T ~ at Mqut, which combined with 
large SSB mass-squared parameters result in b\ and T\ heavier than ~ 3 TeV, much larger than 
the Xi mass. However, dedicated scans can find some solutions where stop co-annihilation does 
occur [47], as shown in Fig. 5a). These solutions tend to be very fine-tuned, since a large value 
of weak scale At must be generated which pushes a normally several- TeV t\ state into the range 
where m,~ ~ m^o . 

b X Al 

We note here that solutions with a neutralino overabundance may be brought into accord with 
the WMAP measured value of CDM by at least two methods. 

1. The Xi may not in fact be the LSP, but instead decays into a (usually) much lighter state. 
This occurs in models with mixed axion/axino (ad) dark matter [61], where Xi ~^ l a - 
Then the neutralino abundance is converted into an axino abundance with [62] ^lah 2 = 
^^Cl^-oh 2 . If rria ~ MeV scale, then the ratio ma/m^o can reduce the putative neutralino 

m -r > Al Al 

mass abundance by many orders of magnitude. In this case of mixed ad CDM, axion- 
domination tends to be favored [63]. 

2. If additional late decaying scalar fields are present in the model, they may get produced at 
large rates via coherent oscillations. If they temporarily dominate the energy density of the 
universe, and then decay to mainly SM particles, they may inject considerable entropy into 
the cosmic soup, thus diluting all relics which are present at the time of decay. Entropy 
injection can occur at large rates for instance from saxion production in the Peccei-Quinn 
augmented MSSM [64, 65], or from moduli production and decay, as is expected in string 
theory [66]. 

In the cases where the neutralino relic abundance is too low - e.g. the low tan (3 higgsino line 
of solutions at ^^»h 2 ~ 10 -3 in Fig. 4 - then the neutralino abundance can be augmented in 
the PQMSSM case where rria > m^o, and additional neutralinos are produced via thermal axino 

Al 

production and decay d — > 7X?, or via saxion cascade decays [65]. In these cases, the CDM tends 
to be neutralino dominated with a small component of axions. 
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Figure 5: Distribution of points from the scan in various planes of sparticles/Higgs masses for 
fj, > and wit = 173.3 GeV. All points satisfy mass bounds, S-physics bounds and have R < 1.05 
and O/i 2 < 0.139. Red and green color represent solutions with tan/3 < 20 and tan/3 > 20, 
respectively. Diagonal black lines represent approximate mass conditions for sizable coannihilation 
or resonance pair annihilation of neutralinos. 
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5 Conclusions 



In this paper, we have used the Isasugra sparticle mass calculator to explore the possibility of b — r 
Yukawa-unified solutions as would be expected to occur in minimal SU(5) SUSY GUT theories. 
Our main assumption is that SU (5) breaks at Q = Mqut to the MSSM as the low energy 
effective theory, so that soft SUSY breaking terms are related by SU(5) boundary conditions at 
the GUT scale. This could occur for instance in 4 — d models with GUT symmetry breaking via 
the Higgs mechanism, or in extra dimensional GUTs where SU(5) is broken via extra-dimensional 
compactification on perhaps an orbifold. We search for sparticle mass spectra which maintain 
fb - fr at M GUT . 

We have found two sets of solutions. At low tan/?, the solutions are characterized by ~ 
f T ~ 0.04, while ft ~ 0.55. These solutions have trig ~ m g ~ 1 — 3 TeV, and only a portion of 
solution space will be accessible to LHC14 searches. A variety of mechanisms are available so 
that thermal production of neutralino CDM may occur in the early universe and generate a dark 
matter abundance at the measured value. However, these low tan f3 Yukawa-unified solutions tend 
to have < 123 GeV, and so may be ruled out if the Atlas/CMS preliminary evidence for a 
Higgs with vrih ~ 125 GeV persists. 

The large tan f3 ~ 35 — 60 solutions exhibit many of the characteristics of t — b — r unified 
solutions which are expected to occur in S'O(IO) SUSY GUTs. Our solutions tend to have fb ~ f T , 
but with a small mis- match with f t . These solutions tend to have multi-TeV matter scalars, but 
nig < 2 TeV. Thus, they relax somewhat the tendency that m g < 500 GeV as occurs in t — b — r 
unified models with gaugino mass unification and fi > 0. This class of models should lead 
ultimately to detection of gluino pair events at either LHC7 or LHC14. They tend to overproduce 
neutralino dark matter except in the (unlikely) case where top-squark co-annihilation may occur. 
The cases of neutralino dark matter overproduction can be brought into accord with astrophysical 
measurements by invoking further neutralino decays (for instance decay to an axino LSP), or via 
entropy dilution of any relics by saxion or moduli decays. 
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